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In heavy-ion fusion reactions involving a well deformed nucleus, one often assumes that the 
orientation of the target nucleus does not change during the reaction. We discuss the accuracy of 
this procedure by analyzing the excitation function of the fusion cross section and the fusion barrier 
distribution in the reactions of ^^''Sm target with various projectiles ranging from ^^C to '"'Ar. It is 
shown that the approximation gradually looses its accuracy with increasing charge product of the 
projectile and target nuclei because of the effects of finite excitation energy of the target nucleus. 
The relevance of such inaccuracy in analyzing the experimental data is also discussed. 

25.60. Pj, 21.60.Ev, 24.10.Eq, 25. 

It is now well established that nuclear intrinsic degrees of freedom strongly influence the fusion cross section in 
heavy- ion reactions at energies near and below the Coulomb barrier Typical examples include the rotational 

excitation of a deformed target nucleus which leads to a large enhancement of the fusion cross section at low energies. 
A characteristic in this case is that the excitation energy of the rotational motion is often much smaller than the 
curvature of the fusion barrier, which determines the time scale of the fusion process. This is the case when, e.g., one 
of the deformed rare earth nuclei or actinides is the target nucleus. One then often calculates the fusion probability 
for each partial wave J following |,| 
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where 9 is the orientation of the deformed target, which will be specified later more precisely. We call this the 
orientation average formula. The upper index S stands for the sudden tunneling approximation. Pj{E,9) is the 
fusion probability for a given orientation. It is determined by solving a one-dimensional Schrodinger equation for 
the relative distance between the projectile and target R, where the potential consists of the nuclear and Coulomb 
components given by 
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with the angle dependent radius of the target nucleus 
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In these equations Vq is the depth parameter of the nuclear potential between the projectile and target, Rp (Rt) 
and Zp (Zt) are the radius and the atomic number of the projectile (target), f3\ is the deformation parameter of the 
target nucleus of multipolarity A, and a the surface diffuseness parameter. The second order term is retained only 
for the quadrupole coupling in the Coulomb interaction. As in almost all analyses of heavy-ion fusion reactions at 
sub-barrier energies, we adopt the no-Coriolis approximation and work in the rotating coordinate frame |^^, where 
the z-axis is taken to be parallel to the radial vector of the relative motion between the projectile and the target 
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nuclei. Also, we assume an axial symmetry for the target nucleus. The parameter 6 is then the angle between the 
symmetry axis and the z-axis in the rotating frame. 

Eq. (|l|) is quite appealing and provides a simple understanding of the large enhancement of the fusion cross section 
at sub-barrier energies. It states that the deformation and the associated rotational excitation of the target nucleus 
result in a distributed fusion barriers, where some of them will be lower than that of the bare potential in the absence 
of these degrees of freedom. The same idea of the orientation average is still correct if one replaces the Gauss integral 
in Eq. (|^) by an appropriate summation even when the rotational band is truncated at a certain angular momentum 
Imax as is always the case in nuclear physics. If only A = 2 is present, the corresponding sum is given by the Imax + 2 
points Gauss quadrature 

The orientation average procedure is exact only in the degenerate spectrum limit, i.e., when the excitation energy 
of the rotational motion is zero, which is not the case in actual nuclei. It is, therefore, important to examine the 
applicability of this formula before one makes a quantitative analysis of the experimental data. A step towards this 
direction has been undertaken in Ref. Q , where an analytic formula has been derived to modify the fusion probability 
in the sudden tunneling approximation by taking the finite excitation energy of the rotational motion into account as 
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where is the excitation energy of the first excited 2+ state of the rotational band. Tq = Tr/il, being the curvature 
of the fusion barrier, is the tunneling time and is a measure of the strength of the channel coupling which causes 
the rotational excitation of the target nucleus during the collision. In Ref. Q|, it was identified with the coupling 
strength at the fusion barrier Rb- 

Eq. (H) is useful to qualitatively discuss the conditions to apply the orientation average formula. It clearly shows 
that the excitation energy of the rotational motion should be much smaller than the curvature of the fusion barrier. 
It also shows that the coupling strength also governs the validity of the formula. Despite these advantages, as we 
see later, Eq. is not tolerable for quantitative discussions. This is partly because Eq. has been derived in 
a perturbation theory, and also by ignoring the radial dependence of the coupling form factor. In heavy-ion fusion 
reactions, the latter is a very crude approximation, because the nuclear and Coulomb couplings interfere leading to 
a strong radial dependence of the coupling form factor Nevertheless, for future reference, we show in Table I the 
exponential factor in Eq. (^, which we call the dissipation factor, for several systems to be examined later. Table II 
lists the bare nuclear potential parameters we used and the Coulomb barrier properties, where the radius parameter 
ro has been introduced as Rp + Rt = ro(Ap ^ + AI/^). 

In this paper we compare the results of the orientation average formula with the numerical solutions of the corre- 
sponding coupled channels equations, which are obtained using the computer code GCFULL joj by keeping the finite 
excitation energy, i.e. 0.082 MeV, of the target nucleus ^^'^Sm. We call the latter as the exact coupled channels 
calculations. We keep only the quadrupole deformation in the target, i.e., A = 2, for simplicity. All the projectiles are 
treated as inert. We truncate the rotational band of ^'''^Sm at lmax='^0'^ member for all reactions. According to the 
table of isotopes the highest member which has been so far experimentally observed is 16+. We have confirmed 
that the comparison between the exact coupled-channels and the orientation average calculations does not almost 
change beyond /„ia2;=12+. We also remark that the coupled-channels calculations almost converge at /mo3;=12+ 
member [|ll| for the and ^^'^Sm reactions, but higher levels introduce non-negligible effects for heavy projectiles. 

Both the orientation average and the exact coupled channels calculations give the fusion cross section which is 
a monotonically increasing function of energy. In order to facilitate to see the accuracy of the orientation average 
formula, we plot in Fig. 1 the ratio of the fusion cross section calculated by solving the exact coupled channels 
equations to that obtained by the orientation average formula. The figure clearly shows that the deviation of the 
results of the orientation average formula from those of the exact calculations gradually increases with the charge 
product of the projectile and target nuclei. This behaviour is consistent with Table I, which has been obtained based 
on Eq. (^, though the actual deviation is much larger than what Eq. predicts. We also note that the deviation 
is significant even for light projectiles at energies below the Coulomb barrier. 

An important question is whether this significantly affects the understanding of the mechanism of heavy-ion fusion 
reactions and the structural informations such as the deformation parameters extracted from the data analyses. In 
this connection, we compare in Fig. 2 the experimental data of the excitation function of the fusion cross section 
and the fusion barrier distribution ||l^,|l^ with the theoretical results calculated by the orientation average formula 
(the dashed line) and by the exact treatment of the corresponding coupled-channels equations (the solid line) for 
^^O-f'^^^'Sm reaction. The second derivative of the cross section times the bombarding energy has been calculated 
with the point difference method with the interval AE = 2 MeV. Similar comparison is done in Fig. 3 for ''*'Ar+^^''Sm 
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reaction. The bare Coulomb barrier is at Vb = 59.41 MeV and 127.57 MeV for ^^O+^^^'Sm and ^''Ar+-'^^'*Sm reactions, 
respectively. In these calculations the deformation parameter P2 of ^^''Sm was determined to be 0.32 by fitting the 
data of ^^'0+^'''^Sm fusion reactions. The same deformation parameter has been used for ^"Ar+^''^"'Sm reactions as 
well, though the effective optimum values can differ in two reactions. Fig. 2 shows that the difference between two 
theoretical lines is much smaller than their deviation from the experimental data. This indicates that one can safely 
use the orientation average formula to data analyses of -"^^O+^'^^^Sm fusion reactions, even though its deviation from 
the exact calculation can be noticeable at low energies as shown in Fig. 1. The top panel of Fig. 3 shows that the 
deviation between the exact and the orientation average calculations is not so drastic even for ^"^Ar+^^^Sm reaction in 
this scmilogarithmic plot. On the other hand, the situation is different for the fusion barrier distribution shown in the 
lower panel. The difference between the two calculations can be more easily recognizable than the case for ^^0+^^''Sm 
fusion reactions. The present theoretical calculations still underestimate the fusion cross sections compared with the 
experimental data for the ''°Ar+^^^Sm reaction at low energies. This will be partly because we ignored the projectile 
excitations, whose excitation energy is as small as 1.46 MeV. 

In summary, we have studied the accuracy of the orientation average formula for fusion cross sections between a 
spherical projectile and a well deformed target by comparing its results with those of the exact coupled channels 
calculations. We found that the results of the orientation average formula significantly differ from those of the 
exact numerical solution of the corresponding coupled channels equations for systems with a large charge product 
of the projectile and target nuclei. This suggests the necessity of the proper coupled-channels calculations beyond 
the orientation average formula for these systems in order to properly identify the role of various channel coupling 
effects and to extract reliable informations on nuclear structure, especially through the analysis of the fusion barrier 
distribution of high precision data. Unexpectedly, we observed a significant deviation of the orientation average 
formula from the exact calculations even for light projectiles in the energy region well below the Coulomb barrier. 
This deviation is, however, much smaller than the deviation of these two calculations from the experimental data. In 
this sense, the orientation average formula is safely applicable to light systems. 
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TABLE I : Dependence of the dissipation factor on the projectile in the fusion of ^^*Sm target estimated based on 
Eq. (5). 



Reactions 


i2c 








^«Si 






4UAr 




0.810 


0.689 


0.567 


0.456 


0.357 


0.274 


0.207 


0.164 



TABLE II : The bare nuclear potential parameters for different systems. 



Systems 


Vo (MeV) 


To (fm) 


a (fm) 


Vb (McV) 


Rb (fm) 


hn (McV) 


12C+^^''Sm 


150.0 


0.950 


1.05 


44.60 


10.81 


3.47 




165.0 


0.950 


1.05 


59.41 


10.81 


3.48 


^"Nc+^-^^Sm 


190.0 


0.945 


1.05 


73.85 


10.87 


3.48 


^^Mg+^'^'^Sm 


225.0 


0.935 


1.05 


88.04 


10.95 


3.50 


^«Si+^''''Sm 


255.0 


0.935 


1.05 


101.58 


11.09 


3.50 


^^S+^^^Sm 


285.0 


0.935 


1.05 


114.89 


11.22 


3.51 


^o^r+i^iSm 


315.0 


0.935 


1.05 


128.00 


11.34 


3.51 


^^Ar+^'^^Sm 


294.0 


0.935 


1.05 


127.57 


11.38 


3.34 
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Figure Captions 



FIG. 1 

The ratio of the fusion cross sections obtained in the exact coupled channels calculations to those obtained in the 
orientation average formula. 

FIG. 2 

Comparison of the experimental data and the results of the orientation average (the dashed line) and the exact coupled 
channels (the solid line) calculations for the ^^0+^^^Sm fusion reaction. The top and bottom panels are the fusion 
excitation function and the fusion barrier distribution, respectively. The data are taken from Ref. |13|| . 

FIG. 3 

The same as Fig. 2, but for the ''*'Ar+^^"'Sm fusion reaction. The data are taken from Ref. [l^ . 
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